Spin Density Functional Based Search for 
Half-Metallic Antiferromagnets 
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We present results based on local spin density calculations of a computational search for half- 
metallic (HM) antiferromagnetic (AFM) materials within the class of double perovskite structure 
oxides LaM M O3 that incorporate open shell 3d (or 4d) transition metal ions M ,M . The pairs 
M'M" = MnCo, CrFe, CrRu, CrNi, MnV, and VCu are studied. La2VMn06 is the most promising 
candidate, with the HM AFM phase more stable by 0.17 eV/cell than the ferromagnetic phase. 
La2VCu06 is another promising possibility, but due to S=i ions quantum fluctuations my play 
an important role in determining the ground state magnetic and electronic structure. This study 
indicates that HM AFM materials should not be prohibitively difficult to find. 



I. INTRODUCTION 

A nonmagnetic metal whose electrical conduction is 
100% spin polarized seems like an oxymoron. However, 
this is only one of the peculiar properties ||l|-^ of a half- 
metallic (HM) antiferromagnet (AFM). A HM AFM is 
also a metal with zero spin susceptibility, a property nor- 
mally associated with insulators, whereas conventional 
metals have x N(E^), the Fermi level density of states. 
Due to mixing of atomic orbitals with neighboring non- 
magnetic ions such as oxygen, magnetic ions commonly 
take on non-integral values of magnetic moment. Yet in 
a HM AFM distinct ions have antialigned moments that 
cancel exactly. |^ The prospect of getting a fully mag- 
netized current from a metal that has no magnetic field 
provides not only fertile ground for research but conceiv- 
able new "spin electronics" devices that rely on the spin 
polarization of the carriers. Recently the likelihood of a 
novel form of superconductivity in HM AFMs has been 
proposed, 0] which further intensifies theoretical interest 
in these unique systems. To date, there is no confirmed 
example of a HM AFM. Is there a real likelihood of dis- 
covering, or even predicting, new HM AFMs, or are they 
destined to remain a theoretical curiosity? 

Characterization of a material as "half-metallic" spec- 
ifies that one spin channel is metallic while the other 
channel is insulating |^ . For a stoichiometric compound 
this results in a spin magnetization that is an integer 
number |0] of Bohr magnetons (/is) per cell. The 
present objective is to predict compounds where M =0 
that have the special properties mentioned above. For an 
initial study one should consider only the simplest case 
of two magnetic ions whose spins (S) are antiparallel. In 
this case the moments will be distinct in shape and ex- 
tent (different spin densities) but will cancel precisely in 
each cell due to the half-metallic nature of the system. 

The only suggestion at present for a HM AFM 
material is the Heusler-like quintinary ordered alloy 



VrMnFesSbyAs. Q Due to the complexity and intricacy 
of the unit cell, and to the intermetallic nature of the 
constituents, it is unlikely to be made in stoichiometric 
form. The perovskite crystal structure AMO3, due to its 
simple crystal structure, potentially very large number 
of members, and strong coupling between magnetic or- 
dering and electronic properties, appears to be an ideal 
system for a search for HM AFM members. We report 
here results of a computational search for candidate HM 
AFMs, based on a linearized augmented plane wave im- 
plementation of spin density functional theory , in the 
double perovskite crystal structure La2M M Oq. More 
than 300 compounds in this structure, shown in Fig. 1, 
have been reported however, very few of these con- 
tain two magnetic ions. A=ha, was chosen for this study 
because of experience with this cation Q and also due to 
the fact that LaMOs compounds in the perovskite struc- 
ture are known to exist for all ions M in the 3d series. 



II. BASIC CONSIDERATIONS OF THE SEARCH 

Magnetic 3d ions are characterized by a nominal d" 
configuration where 1 < n < 9. Since each such config- 
uration can be realized by more than one ion (in differ- 
ent charge states), and there are several uni-, di-, and 
trivalent cations A to choose from, there are thousands 
of magnetic double perovskite compounds that might 
be considered. Although attributes of a material that 
make it a better candidate comprise a substantial list, 0] 
there is one overriding requirement: the moments must 
be equal in magnitude so they may cancel to give A^=0. 
This requirement narrows the number of choices consid- 
erably. 
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FIG. 1. The double perovskite crystal structure. The black 
and white spheres are transition metal (M', M") ions, the 
gray oxygen ions form an octahedron around each metal ion, 
and the cation (viz. La, but not shown), lie between eight 
MOe octahedra. The figure illustrates that the octahedra 
around different M ions are allowed by symmetry to be dif- 
ferent sizes. 

A. Combinatorics of A'A"M'M"Oe, 

To illustrate the combinatorics in the class of com- 
pounds A2M'M"0e, we consider the number of double 
perovskite structure compounds that might be relevant 
for our search. There are nine possible open configu- 
rations, giving distinct pairs of 3d ion. However, id 
ions are also likely to be magnetic in these compounds, 
and including them the number of possibilities becomes 
pairs. Since an ion can be in a few different charge 
states, each 3d" or 4d" configuration can be attained by 
more than one ion; we take as an average two charge 
states per ion, giving another factor 2^. The cation A 
can be chosen from di- and tripositive cations (including 
the rare earths) and even some univalent ions, amounting 
to some 25 ions. The number of compounds then is of 
the order of 

25 X 2^ X 15000. (1) 

2 ^ ' 

Considering also the possibility of splitting the cations 
A — » A' A" leads to an addition factor of or a total of 
the order of 2x10^ possibilities. 

This vast number reduces drastically if one considers 
the requirements imposed by the HM AFM state. Most 
crucially, the magnetic moments of the ions must be equal 



in magnitude, so they have the possibility of cancelling. 
To do this, one separates the ions into the five classes 
S=^, 1, |, 2, |, and only pairs with equal spins need be 
considered. Supposing the ions separate evenly into these 
classes (not really the case) with about four per class, 
the factor in Eq. (1) drops to 5 x The cation 

charge further restricts the pairs of ions. Let us consider 
only the case A=La!^^ , which is all we consider in this pa- 
per anyway (removing the factor of 25 in Eq. (1)). Spec- 
ification of this trivalent cation implies that the charge 
states of the two d cations must sum to six, effectively 
limiting the choice to (3+,3+) or (2+,4+) pairs ((l-|-,5-|-) 
is much more rare). Due to complexities in trying to pre- 
dict the spins of each ion (which may be environment de- 
pendent, see below) we will not attempt to enumerate the 
possibilities, but the number of possible La-based double 
perovskites is possibly no more than 25. 

B. Estimation of Stable Moments 

To use the constraint of equal moments it is necessary 
to be able to make a realistic prediction of the spin mo- 
ment of a given ion in the crystal. For a magnetic ion 
with an open 3d shell, the moment m (in /is) is 

m = snjj^s = n+—n-, (2) 

a.s 

where s=±l denotes the spin direction, ria^s is the occu- 
pation number of 3d crystal field level a = or t2g , with 
spin s, and rig is the total s occupation. The energy of 
an ion in the perovskite octahedral field is 

^ / 3 s 

Eron^ '^\ne^,s{-^^CF - ^Istm) (3) 
s 

2 s \ 

+ nt^^^^s{-^^CF - ^Istm)y (4) 

where Ac_f is the crystal field splitting between the Cg 
and t2g levels, 1st is a Stoner-like parameter giving the 
exchange (magnetic) energy ±iAea: 

Aej; = Ist-m (5) 

proportional to the moment m on the ion, and n = + 
rt_. This energy can be expressed as 

Eion = Ac_F ^ Tis i/s - i/stm^, (6) 

s 

Ds=^ na,s Val ^Ha.s. (7) 

Q a 

Here v y& f , — | for eg,t2g, respectively, and this formula 
for the energy illustrates the balance between crystal field 
energy and magnetic energy that must be minimized to 
obtain the stable configuration of the ion. 
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The ratio 7 = Igt/^CF determines whether the ion 
takes on a high spin or low spin configuration. S=| is 
attained only by a cf^ ion (viz. Mn^+, Fc'^+, or Co +). 
S=2 occurs only for df^ and ions, and S=| only for 
and ions. S=l is attained not only by d^ and d^ high 
spin ions, but also by the low spin d'* ion (i2g,+*23,-) if 7 
is small. S=^ occurs for low spin d^ and d"^ ions as well 
as for the single electron d^ and single hole d^ ions. 

These guidelines hold in the strong crystal field limit. 
Within a crystal there are additional effects. The ki- 
netic energy, which results from electron hopping from 
ion to ion {via the intervening oxygen ion in the per- 
ovskite structure) drives the moment away from integer 
values and has very substantial effects in perovskite ma- 
terials, structural distortions that lower symmetry and 
alter hopping amplitudes. There is in addition the pos- 
sibility of ions changing their charge state (determined 
by differences in site energies and by intra-atomic repul- 
sion). However, these single ion energies, with 1st and 
AcF given from density functional calculations or from 
experiment, have been used to guide our initial choice 
of ion pairs that we have examined more closely using 
self-consistent calculations. 

Fig. 2 indicates the behavior of the ion energy in the 
strong crystal limit for two values of 7, for each integer 
occupation n as the moment m varies over allowed values 
[m < n for n <5, m < 10 — 71 for 5< n <10). All oc- 
cupations Ua^s were varied within the limits of their con- 
straints until the energy was minimized for fixed n and m 
(for the isolated ion, and neglecting spin-orbit coupling, 
m=2S). The energy reflects the downward parabolic en- 
ergy gain from polarization, together with a linear in- 
crease in energy when one crystal field level is filled and 
a higher one begins to become occupied. There is a 
strong tendency toward the "high spin" (largest possible 
m, Hund's rule) configuration; however, for large Ac_f 
relative to 1st a low spin configuration can become sta- 
bilized: note the n=6 curves in both plots, where a non- 
magnetic 771=0 result is stable or metastable. Low spin 
configurations may also occur for 7i=4, 5, or 7, as noted 
above. The curves of Fig. 2 are only guidelines; fully self- 
consistent local spin density calculations reported below 
vary both 71 and 711 for a given spin alignment until the 
energy in minimized. For the perovskite structure, values 
of AcF are 2-3 eV. 



III. RESULTS FOR THE CHOSEN COMPOUNDS 

The double perovskite compounds that have been 
studied ||^ are listed in Table I. We discuss them in turn. 
Note that, while we denote global spin directions by up 
(-|-) and down (-), the terms 'majority' and 'minority' 
refers to each specific ion (more occupied, or less occu- 
pied) and does not specify a direction of the spin. Also, 
due to strong hybridization with the oxygen 2p states, it 




Magnetic Moment (/x ) 




Magnetic Moment (/x^) 
FIG. 2. Energy of a magnetic ion in a crystal field Acf, 
versus magnetic moment, for various ionic configurations d". 
(a) 7=0.1. (b) 7=0.2; the 'spin only' dotted curve indicates 
the asymptote as the crystal field term vanishes.. Each curve 
is the sum of a piecewise linear increasing crystal field term 
and a negative quadratic magnetic term. 



is not possible even to define an ionic moment precisely. 
Values given in Table I and below provide the nominal 
charge and an estimate of the effective ionic moment, 
which often is not very near any integral value. 



A. Mn^+,Co 



d\d'' 



These ions were chosen as a likely S=2 pair. In the 
ionic picture, the Mn'^"'" ion is expected to be high spin 
(4 ^b), while Co'^"'" may be high spin (4 /is) or non- 
magnetic, see Fig. 2. In fact, nominally trivalent Co ions 
are known sometimes to occur in a low but non-zero spin 
(1-2 /Is) state as well. When the Mn and Co moments 
are antialigned, indeed a HM AFM state is obtained, al- 
though the moments are only « 2.8/iB. The densities 
of states (DOS), shown in Fig. 3(a), indicate that the 
ionic picture is followed closely (Ac_f=1.5-2 eV, Aga, w2 
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TABLE I. Local spin density functional results for mag- 
netic double perovskites La,2M' M"0(i. Approximate cal- 
culated d occupations n are given as d". Approximate 
spin-only moments m(M'), m(M") (in ps) are not always 
near an integer value, due to strong hybridization with oxy- 
gen. Half-metallic character is denoted by an integer value 
of Mtot, and HM AFMs are denoted by AFM^ Relative 
energies are per transition metal ion, with the first quoted 
magnetic alignment taken as the reference. FiM = ferrimag- 
netic. 



M' 


M 


Order 


m(M ) 


m(M \ 


Mtot 


Ere;(eV/M 


Mn^+ d^ 


Co^+ d' 


AFMf 


2.8 


-2.8 


0. 


0.00 






FM 


3.3 


1.3 


4.60 


-0.46 




Fe^+ d^ 


FiMl 


2.6 


-0.6 


2. 


0.00 






FiM2 


2.0 


-4.0 


-2. 


-0.12 






FM 


3.0 


4.0 


7.15 


-0.03 


Cr'+ d' 


Rn'+ d^ 


FiM 


2.5 


-0.5 


2. 


0.00 






PM 








+0.Q7 




Ni^+ d« 


FiM 


2.0 


-1.4 


0.60 


0.00 






FM 


2.4 


1.6 


4. 


-0.15 


Mn'+ d^ 


d^ 


AFMl" 


1.6 


-1.6 


0. 


0.0 






FM 


1.9 


1.1 


3.00 


-1-0.17 


d^ 


Cn^+ d' 


AFMl" 


0.7 


-0.7 


0. 


0.00(!) 






FM 


0.7 


0.7 


1.38 


0.00(!) 






PM 








-H0.06 



eV). The conducting channel has 75% Co d, 25% Mn d 
character at the Fermi level. When the spins are aligned, 
a high-spin Fe and low-spin Co result was obtained, re- 
flecting the strong difference in hybridization that arises 
due to the type of magnetic alignment. The spin-ahgned 
phase is 0.46 eV per M ion lower in energy than the HM 
AFM phase, so the sought-after HM AFM phase is at 
best metastable. 



> 

73 
CD 

o 

0) 

P 



10 

5 



-10 
10 



10 
10 



5 

-5 
-10 



La^VCuC 


1 

1 


+ spin 

a: 


(a) 


Cu:e /-'' V:e 
Ig - spin 


1 1 1 1 , 1 1 1 


La MnVO„ 


-1- spin 

Mn:t3^ 


Mn:t^ 

(b) 


1 Mn:e ' ag 
E ' - spin 

r 




La,MnCoO„ 
Mn:t,^ 


'■■ Co:t + spin 

',1 2s 

'f; Mn:e 


(c) 


Core 1 

' 1 


YMn:t,^ . 
1 — spin 



Energy (eV) 

FIG. 3. Atom-projected densities of states for the three 
HM AFM states found in the double perovskites: (a) Mn-Co; 
(b) Mn-V; (c) V-Cu. Note that in (a) and (c) the t2g and eg 
states are distinguishable for each ion and each spin, whereas 
this is not the case for the strongly mixed metallic channel in 
(b). 



B. Cr^+,Fe^+: d^d^ 

This pair was chosen in the anticipation that differing 
charge states [Cr^+,Fe^+ : d*,d'^] might result, leading 
to either an S=2 pair or an S=l pair. However, charge 
differentiation did not occur. Two distinct solutions for 
antialigned moments have been obtained, one has high 
spin Fe (4 fis) parallel to the net moment, and the other 
has low spin Fe ("1 /^b") antialigned with the net mo- 
ment. Both are HM; however, they are ferrimagnetic 
(FiM) rather than AFM. A FM, high spin solution also 
was found. Of these three spin configurations, the high Fe 
spin FiM state is lowest in energy. All three states can be 
characterized as Cr'^+,Fe'^+ (d'^,d^) whose moments will 
not cancel. Thus current evidence is that this system is 
not a prospect for HM AFM. However, a HM FiM state 
is hkely, and is itself worthy of study. 



C. Cr3+,Ru^+: d^ d^ 

The Ad Ru ion is isoelcctronic with Fc, so this pair pro- 
vides another possibility for a S=2 or S—1 pair with un- 
equal charges (Cr^+,Ru'*+). The unexpected behavior in 
Ru-based perovskite compounds (ferromagnetism (FM) 
in SrRuOs @J§1; superconductivity at IK in Sr2Ru04 
Q) also prompted us to check the pair Cr-Ru. How- 
ever, like Cr-Fe, the Cr-Ru pair are predicted to be iso- 
valent, and Ru is found to be magnetic with a low spin. 
With antialigned moments, the compound is a HM FiM 
with a net moment of 2 fis- Using fixed spin moment 
methods it is possible to force the net moment to vanish 
slowly during continued self-consistency. Sometimes this 
results in the discovery of another phase, usually one that 
is metastable. In this case, however, forcing the net mo- 
ment to vanish resulted in the destruction of the moment 
on both ions, at the energy cost of 0.67 cV/ion, rather 
than resulting in a (possibly HM) AFM state. As in all 
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of these compounds, this pair may have magnetic states 
that we have not located in these calculations. 

D. Cr=^+,Ni=^+: d^ 

This pair was chosen as a likely S=| case. For small 
crystal field, the dJ ion may be high spin (3 hb), thereby 
balancing the cfi moment. For these pairs the antialigned 
configuration resulted in a FiM net moment of 0.60 /is, 
near but not at the HM AFM result. The FM alignment 
is lower in energy by 0.15 eV per ion. This is a case where 
relaxation of the structure (volume and/or distortion) 
might result in a HM AFM phase. 

E. Mn='+,V='+: d'^,d'^. 

This choice, if Mn^+ assumes low spin, is a possibility 
for an S=l pair. A low spin indeed resulted, for both 
antialigned and aligned moments. The DOSs are shown 
in Fig. 3(b); we discuss the unusual structure of the 
metallic channel below. The antialigned ordering results 
in a HM AFM state, moreover this state is 0.17 eV/ion 
lower in energy than the FM alignment. Thus this pair 
provides a strong candidate for the HM AFM that this 
search hoped to locate. This compound is discussed fur- 
ther in the next section. 

F. V*+,Cu2+: d^,d\ 

This was the choice for an pair. This pair is 

unique in our studies to date, in that differing charge 
states actually are obtained. The differing ionic radii 
as well as the differing charge states of this pair of ions 
should give a strong preference for well ordered struc- 
tures. Both aligned and antialigned moment solutions 
were obtained, with identical energies. Moreover, the an- 
tialigned phase is a HM AFM, whose DOS is compared 
in Fig. 3(c) with the other HM AFM states. Due to the 
identical energies of the different alignments and their 
small spins, this system appears to be a strong candi- 
date for a three dimensional quantum magnetic system 
as T^O {viz. quantum spin liquids, heavy fermion met- 
als, Kondo insulators), rather than a simple HM AFM or 
an exotic superconductor. 

G. Synopsis 

Of these six pairs of transition metal ions giving the 
double perovskite structure, three have led to at least 
a metastable HM AFM phase. This amount of success 
is remarkable considering there was previously no viable 
candidate. One of these is clearly only metastable (the 



Mn-Co pair) while Mn-V and V-Cu in the LaM' m" O3 
compound are not unstable towards ferromagnetic align- 
ment of the moments. 

The calculated DOSs of these three HM AFM states, 
shown in Fig. 3, reveal that qualitatively different type 
of gaps can occur in the insulating channel. Correlated 
insulating oxide compounds are identified as Mott insu- 
lators, if the gap lies between d states on the metal ion, 
or as charge-transfer insulators if the gap lies between 
occupied O p states and metal ion d states With- 
out experimental input we cannot judge the strength of 
correlation effects in these compounds. However, adapt- 
ing this terminology to the present materials, the Mn-V 
compound (Fig. 3(b)), with a gap between Mn t2g and 
Cg states, is the analog of a Mott insulator. The Mn- 
Co and V-Cu compounds, on the other hand, are inter- 
transition-metal charge transfer insulators, with the gap 
lying between occupied d states on one metal ion and 
unoccupied d states on the other metal ion. 

IV. DISCUSSION 

Although the local spin density functional calculations 
used here generally give good charge densities and in a 
majority of systems (including perovskites Q) predict 
magnetic properties as well, several questions remain. 
Whether the proposed compounds can be made can be 
answered only by experiment (competing phases are too 
numerous to calculate). Magnetic moments, type of spin 
ordering, and ionic charges can be sensitive to volume 
and most calculations reported here were carried out only 
at the representative cubic perovskite lattice constant of 
3.89 A. However, variations of the volume for the Mn-V 
and V-Cu compounds confirmed that this volume is re- 
alistic (near the minimum of energy) and that the HM 
AFM phases persisted at nearby volumes. Relaxation of 
the positions of the O ions as allowed by symmetry (see 
Fig. 1) was carried out for the Mn-V compound. The 
oxygen octahedron relaxed inward around the smaller Mn 
ion ||ll[] by only 0.02 A, and the frequency of O vibra- 
tion about this minimum (an observable Raman-active 
Ag eigenmode) is calculated to be 465 cm~^. This relax- 
ation actually stabilizes the HM AFM state; conversely, 
displacing the oxygen octahedron toward the V ion (by 
0.08 A) drives the compound from the HM AFM phase 
to a FiM phase. 

Possible correlation effects deserve more comment. A 
limitation of local spin density functional calculations 
used here is the inability to predict whether transition 
metal compounds such as these are correlated electron 
systems. In the La2VCu06 compound, for example, the 
bandwidth of the conduction band of the HM AFM state 
is only 1 eV wide, and strong on-site repulsion ("Hubbard 
U") of electrons will tend to drive the metallic channel 
insulating. If this occurs, it may provide an example of 
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FIG. 4. Isosurfaces (at ±0.003 a.u.) of the spin density 
in the HM AFM state of La2VCu06, illustrating the very 
different spin densities on the two magnetic ions. Two double 
perovskite cells are shown; the complete structure is cubic. 
Upper: the largest (light shaded) surface surrounds the V ion. 
The O ions, which appear as dumbbells (truncated by the cell 
edges in some cases) are polarized in the same direction as the 
smaller (dark) Cu ion. 

yet another new phenomenon in these systems: a case 
in which one spin channel is a Mott insulator while the 
other is better described as a band insulator, or possibly a 
generalized charge-transfer insulator. An extension of the 
local spin density method, referred to as the LSDA+U 
method, often provides an improved mean-field de- 
scription of correlated insulators. However, at present it 
is inapplicable to metals so its use is restricted to mate- 
rials that are known to be insulators. 

In La2MnV06, on the other hand, the bandwidth of 
the metallic channel is more than 3 eV wide, and corre- 
lation effects should be much less important. The dis- 
tinctions in Fig. 3(b) are easy to understand. The Mn 
majority t2g states lie lower than any V states, and also 
hybridize less weakly with O 2p states than do Cg states. 
Hopping within this band therefore must go Mn-0-0- 
Mn, hence the band is quite narrow (1 eV). The minority 
Mn t2g states lie in the same range as the majority V t2g 
states; however, since the ionic moments are antiparallel 
these states have the same spin direction (-1- in Fig. 3(b)) 
and can form a mutual, relatively broad, band based on 
Mn-O-V-O-Mn-... hopping. Occupation of this compar- 
itively broad band results in cancellation of part of the 
Mn moment, which by Eq. (3) results (self-consistently) 
in reduced exchange splitting and leads to the low spin 
Mn ion. Hence it is reasonable to expect the conducting 
character of this band to survive correlation corrections. 

An example of a spin density for a HM AFM phase 



is shown in Fig. 4 for the V-Cu compound. The spin- 
density isosurfaces illustrate very graphically the differ- 
ence between the up and down spin densities. In fact, un- 
like conventional perovskite antiferromagnets where the 
oxygen is polarized only in a dipole form with no net mo- 
ment, here the O ion has a net moment that lies in the 
same direction as the Cu spin. This unusual form of spin 
density for an AFM should be more readily apparent in 
the magnetic form factor measured in polarized neutron 
scattering experiments than in the typical case in which 
there is no net moment on the O ion. The O moment in 
this case is only of the order of 0.01 /is, however. 

This computational search indicates that the double 
perovskite class of compounds provides a fertile environ- 
ment for half-metallic antiferromagnets, a new type of 
magnetic material for which there were previously no 
viable possibilities. Although this search was confined 
to the La cation (which may be considered representa- 
tive of trivalent cations), mixed cation compounds such 
as A^'^ B^^ M' M"Oq show strong tendencies to form or- 
dered structures Q and may also provide good candi- 
dates. This work should not be interpreted as suggesting 
that perovskites provide the most likely possibility; cer- 
tainly other crystal structures deserve strong considera- 
tion. The point here is that it should not be a formidable 
problem to fabricate HM AFM compounds, so the study 
of their myriad unusual properties may commence. 
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